Fission half-lives of super-heavy nuclei in a microscopic approach 
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A systematic study of 160 heavy and super-heavy nuclei is performed in the Hartree-Fock- 
Bogoliubov approach with the finite range and density dependent Gogny force with the DIS pa- 
rameter set. We show calculations in several approximations: with axially symmetric and reflexion 
symmetric wave functions, with axially symmetric and non-reflexion symmetric wave functions and 
finally some representative examples with triaxial wave functions are also discussed. 

Relevant properties of the ground state and along the fission path are thoroughly analyzed. 
Fission barriers, Qa-factors and lifetimes with respect to fission and a-decay as well as other 
observables are discussed. Larger configuration spaces and more general HFB wave functions as 
compared to previous studies provide a very good agreement with the experimental data. 



PACS numbers: 21.60.Jz, 21.10.Dr, 21.10.-k, 21.10.Pc, 25.85.Ca, 27.90.-fb 



I. INTRODUCTION 



The stability and structure of nuclei at the upper end 
of the nuclear chart is a hot topic in contemporary nu- 
clear physics. It is a challenging task to answer which 
heavy nuclides may exist and what properties they may 
have. Therefore strong efforts have been made in the 
experimental developments as well as in the theoretical 
description of super- heavy elements (SHE's). 

In the last decades a huge progress in the synthesis 
of new elements has been achieved in world leading lab- 
oratories like the GSI, Darmstadt [1-8], JINR, Dubna 
[9-14], and RIKEN, Tokyo [15-18]. The first SHE's with 
Z < 113 and iV < 165 were produced in cold fusion reac- 
tions. In these experiments, involving neutron rich pro- 
jectiles and spherical targets (■^°^Pb or ^"^Bi), weakly ex- 
cited compound nuclei were produced which cooled down 
by the emission of only one or two neutrons. Further 
experimental progress in the synthesis of the heaviest el- 
ements was achieved by hot fusion reactions in which 
targets of deformed actinide were bombarded with the 
doubly-magic nucleus ^^Ca. The compound nucleus cre- 
ated in this way was more excited and three or more 
neutrons were emitted. These reactions succeeded in the 
synthesis of new elements up to Z = 118 and = 176 
[10]. The first observation of the element 117 was possible 
lately [13] using a radioactive ^^^Bk target. New possi- 
bilities for the synthesis of new isotopes will be opened in 
heavy ion collisions with radioactive ion beams [19, 20]. 
Nowadays, many other laboratories are involved in the 
exploration of SHE's. Thus the experimental groups from 
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Berkeley [21], GANIL [19], Livermore [11, 22], Jyvaskyla 
[23, 24], and Oak Ridge [14] are working in this direction. 
They would bring in the nearest future further informa- 
tion on the stability and properties of the SHE's and an 
independent verification of the existing data. 

In a parallel way to the experimental efforts, the prop- 
erties of the she's have been also investigated in various 
nuclear models. A proper description of trans-fermium 
nuclei is a great challenge for any theoretical model. Usu- 
ally, the parameters of the theories on atomic nuclei are 
adjusted to the stable isotopes and then extrapolated to 
the region of heavier systems. Therefore many tries in 
different theoretical approaches are performed to fore- 
see the stability and the structure of the heaviest nuclei. 
A detailed review of the theoretical analysis of SHE's 
can be found in Ref. [25]. The first theoretical inves- 
tigations on the stability of heavy nuclei were made in 
the 60's. It was noticed that shell effects could stabilize 
nuclei heavier than those known at that time [26, 27]. 
Calculations made in the macroscopic-microscopic model 
with Strutinsky shell correction predicted the "island of 
stability" . Large values of shell energies were obtained 
a,t Z = 108, N ~ 162 for prolate deformed nuclei and 
Z = 114, iV = 184 for spherical ones [26, 28]. In the 
last decades many calculations have been made provid- 
ing more and more precise predictions. The fission bar- 
riers and the ground state properties were calculated us- 
ing macroscopic-microscopic methods with large range 
of deformation parameters and nuclear shapes, including 
reflection and axial symmetry breaking [29-40] . 

Self-consistent methods also provided many results on 
flssion barriers and half-lives. Important results have 
been obtained in the relativistic mean-field (RMF) [40- 
43] , the Hartree-Fock (HE) approach with Skyrme forces 
[41, 42, 44-48] and the Hartree Fock-Bogoliubov (HFB) 
theory with Gogny forces [49-51]. The first calculations 
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of fission barriers were performed in the axial and reflec- 
tion symmetric regime but later all relevant deformations 
were considered in the minimization of the energy. 

It is well known that the liquid drop model docs not 
predict a fission barrier in the heaviest nuclei and that 
the stability of the trans-fermium nuclei is achieved by 
the shell effects. The self-consistent quantum mechani- 
cal methods (RMF, HF and HFB) as microscopic theories 
are the perfect tools for the analysis of SHE's. Moreover, 
in the self-consistent calculations all possible shapes of a 
nucleus are considered in the minimization process. In 
contrast, the most of the macroscopic- microscopic mod- 
els are restricted to some pre-defined classes of deforma- 
tions and only "optimal shapes" [52] allow to obtain any 
configuration of a nucleus. Therefore they are very suit- 
able to describe large deformations of nuclei around the 
scission point. A degree of freedom which plays an im- 
portant role along the fission path are the pairing correla- 
tions [53] . Since in the HFB theory the particle-hole and 
the particle-particle matrix elements are treated on the 
same footing, the proper consideration of pairing along 
the whole fission path is guaranteed. This method has 
been successfully applied in many aspects of low energy 
nuclear physics, in particular in the description of fission 
barriers of heavy nuclei [49, 50, 54] . Another theoretical 
quest was to discover a semi-empirical formula describ- 
ing a emission half-lives [55, 56]. These investigations are 
very important as a radioactivity is the dominant decay 
channel in many SHE's. 

The purpose of this article is to perform a systematic 
study of SHE's with respect to their stability and ground 
state properties in the framework of the HFB theory with 
the density dependent finite range Gogny force and the 
DIS parametrization. In our analysis we include the re- 
gion of the well known Fermium (Fm, Z = 100) and 
Nobelium (No, Z = 102) elements to compare our pre- 
dictions with the available experimental data. We show 
results for the heavier even-even nuclei Rutherfordium 
(Rf, Z = 104), Seaborgium (Sg, Z = 106), Hassium (Hs, 
Z = 108), Darmstadtium (Ds, Z = 110) and Coperni- 
cium (Cn, Z = 112, which was named two years ago 
[57]). Heavier elements with Z = 114 — 124, without 
given name so far, are also considered. We limit our 
study to < 190 isotopes. These nuclei are nowadays 
in the main stream of interest of experiments with SHE's. 

A large amount of information on nuclear structure 
and stability can be obtained from properties of the 
ground states of the SHE's. Consequently, we start our 
investigation with the description of the ground state 
characteristics. Deformations, pairing energies and two- 
nuclcon separation energies are analyzed and collated 
with the single-particle energy level scheme. The ground 
state energies can be used to evaluate the Qa values 
which are necessary to calculate the probability of a emis- 
sion - one of the dominant decay modes in SHE. A com- 
peting process to a decay is the spontaneous fission. To 
analyze this mode we determine the fission barriers for 
all mentioned SHE's as a function of the quadrupole mo- 



ment Q2. The calculations were performed in an axial 
basis, although we are aware of the non-axial effects on 
the height of the barrier and we discuss them in a few 
selected cases. The impact of the octupole deformation 
on the potential energy along the fission path is crucial 
in the determination of fission barriers of SHE's. This is 
taken into account by allowing non-reflexion symmetri- 
cal shapes. The calculations were performed in a large 
deformed harmonic oscillator basis paying special atten- 
tion to the proper optimization of the oscillator lengths 
and to the convergence of the calculations with the size 
of the basis. Next, using the WKB approximation, we 
calculate the fission half-lives. Finally, the comparison of 
the half-lives for a decay and spontaneous fission allows 
to predict the stability of the heaviest nuclei. 

This article is organized as follows. In Sec. II we briefly 
describe the constrained Hartree-Fock-Bogoliubov calcu- 
lations. The description of the ground-state properties of 
SHE's is shown in Sec. Ill, flssion barriers in Sec. IV. and 
half-lives of SHE's are discussed in Sec. V. Finally, Sec. 
VI contains a summary and some concluding remarks. 



II. THEORETICAL MODEL 

In our research we will apply the self consistent 
Hartree-Fock-Bogoliubov theory with the finite range 
density dependent Gogny force. In the numerical ap- 
plications we use the DIS [58, 59] parameterization of 
the Gogny interaction. The DIS parameters were ad- 
justed [58] to give a better surface energy term (crucial 
for a proper description of the fission phenomenon). The 
choice of the Gogny force with the DIS parameteriza- 
tion is based on the fact that whenever this interaction 
has been used to describe low energy nuclear structure 
phenomena an, at least, reasonable agreement with ex- 
periment has been obtained. This degree of agreement 
has been obtained both for calculations at the mean field 
level and beyond [60-76]. 



A. Details of self-consistent HFB calculations 

In the microscopic HFB calculations we have used the 
computer code presented in Ref. [78], see also [54] where 
special attention was paid to an accurate computation 
of the matrix elements of the Gogny interaction for very 
large basis like the one used in this paper. The self- 
consistent equations have been solved by expanding the 
quasiparticle creation and annihilation operators on finite 
bases of axially symmetric deformed harmonic oscillator 
(HO) cigenfunctions. The size of the bases used depends 
upon two parameters, Nq and q, which are related to the 
allowed range of the HO quantum numbers through the 
relation 

-ri;, + (2nj_ -I- jmj) < A^o- 
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Along the perpendicular direction we take iVo shells, 
(i.e. 2nx + |m| = 0, . . . , Nq) and along the z direction 
we include up to qNq shells depending on the value of 
2n_L + |to|. In the present study we have used q = 1.5, a 
value which is suited for the elongated shapes along the 
z direction typical of the fission process, and A^o = 15. 
Other parameters characterizing the HO bases are the 
oscillator lengths b± and bz. These two quantities have 
been determined, for each calculated wave function, as 
to minimize the HFB energy. In order to study the dif- 
ferent paths to fission, in our calculations we have used 
as constraints the axial quadrupole {Q2) and octupole 
{Qs) moments, with Qx = r"^ Px{cos{6)) . Higher mul- 
tipolarities are adjusted in the self-consistent process to 
minimize the energy. To study the impact of triaxiality 
effects we have also carried out calculations for a few nu- 
clei where the axial symmetry requirement was released 
but the left-right symmetry was imposed. 

In the calculations the Coulomb exchange energy has 
been treated in the Slater approximation [80, 81]. Ad- 
ditionally, the Coulomb and the spin-orbit contributions 
to the pairing field have been neglected. Finally, the 
two body kinetic energy correction (2b-KEC) is not in- 
cluded in the variation process because, for heavy nuclei, 
it remains almost constant for most of the physical con- 
figurations. As this term was included in the fitting of 
the force, we have to include its contribution at the end 
of the calculation in order to obtain reasonable binding 
energies. See also [54, 76] for a detailed discussion of the 
relevance of neglected or approximated terms in these 
calculations. 

We have also subtracted from the HFB energy the ro- 
tational energy corrections (REC) stemming from the 
restoration of the rotational symmetry. This correction 
has a considerable influence on the energy landscape (and 
therefore on the height of the fission barriers) as it is 
somehow proportional to the degree of symmetry break- 
ing and therefore proportional to the quadrupole mo- 
ment. A full calculation of the REC would imply an 
angular momentum projection [69, 83, 84] which is only 
feasible for light nuclei. In order to estimate the REC 
we have followed the usual recipe [85] of subtracting to 
the HFB energy the quantity (AJ^)/(2Jy), where {A,P) 
is the fluctuation associated with the angular momentum 
operators in the HFB wave function and Jy is the Yoccoz 
moment of inertia [85] . This moment of inertia has been 
computed using the "cranking" approximation in which 
the full linear response matrix appearing in its expres- 
sion is replaced by the zero order approximation. The 
effect of the "cranking approximation" in the Yoccoz mo- 
ment of inertia was analyzed with the Gogny interaction 
for heavy nuclei in [86] by comparing it with the one ex- 
tracted from an angular momentum projected calculation 
(see also [83] for a comparison in light nuclei). The con- 
clusion is that the exact REC is a factor 0.7 smaller than 
the one computed with the "cranking" approximation to 
the Yoccoz moment of inertia for strongly deformed con- 
figurations (a similar behavior has been observed for the 



Thouless-Valatin moment of inertia in [87]). We have 
taken this phenomenological factor into account in our 
calculation of the REC. 



B. Evaluation of lifetimes 

The evaluation of the spontaneous fission half-life is 
carried out in the standard WKB framework where Tgf 
is given (in seconds) by [88] 



Tsf = 2.86 ■ 10-21(1 + exp(2S')) . 



(1) 



In this expression S is the action along the Q2 con- 
strained path which is given by 

rb 

S= / dQ2V'^B{Q2){V{Q2) - Eo) . (2) 

a corresponds to the Q2 value of the ground state and b 
the Q2 value where the potential energy equals the one 
of the ground state. For the collective quadrupole inertia 
B{Q2) we have used the adiabatic time dependent HFB 
(ATDHFB) expression computed again in the "cranking" 
approximation and given by [89] 



M_3(Q2) 

Batdhfb(Q2) = J^T^^ 



with 



M_„(Q2) = 



{E^ + E^Y' 



(3) 



(4) 



Here is the two-particle-zero-hole component of the 

quadrupole operator Q2 in the quasiparticlc representa- 
tion [85] and E^ are the quasiparticle energies obtained 
in the solution of the HFB equation. 

In the expression for the action the collective potential 
V{Q2) is given by the HFB energy (with the 2b-KEC 
and REC corrections) minus the zero point energy (ZPE) 
correction eo((32) associated with the quadrupole motion. 
This ZPE correction is given by 



eo(<32) - 2'^('52)-Batdhfb(<32) 



where 



G{Q2) 



M_2(Q2) 

2Ml^{Q2) 



(5) 



(6) 



Finally, in the expression for the action an additional pa- 
rameter Eq is introduced. This parametcir can be taken as 
the HFB energy of the (metastablc) groimd state. How- 
ever, it is argued that in a quantal treatment of the prob- 
lem the ground state energy is given by the HFB energy 
plus the zero point energy associated to the collective 
motion. To account for this fact, the usual recipe [90] is 
to add an estimation of the zero point energy to the HFB 



4 



energy in order to obtain Eq. In our calculations we have 
taken a zero point energy of 0.5 MeV for all the isotopes 
considered. 

In some SHE's around N = 176 a weakly oblate de- 
formed ground state can be found. The fission path in 
these nuclei obviously does not go along axial symmet- 
ric shapes through the spherical configuration. Such a 
nucleus rather takes triaxial shapes to reach the pro- 
late saddle point with the similar absolute value of the 
quadrupole moment. As the energies on the triaxial part 
of the fission barrier are very small in comparison with 
the saddle point [91] they will contribute insignificantly 
to the action integral. Therefore we will neglect them in 
the calculation of half-lives. 

To calculate the a-decay half-lives we use the phe- 
nomenological formula of Viola and Seaborg [92] 

logr„[yr] = {aZ + 6)(Q„)-i/2 ^ ^ ^) _ 7 5 (7) 

with Z the atomic number of the parent nucleus. The 
Q-factor of the decay, Qa, is obtained from the calcu- 
lated ground state binding energies with the experimen- 
tal value E{2, 2) = -28.295674 MeV [93] 

Q^iZ, N) = E{Z, N) -E{Z-2,N-2)- E{2, 2) (8) 

The constants in Eq. (7) are: a = 1.66175,6 = 
-8.5166, c = -0.20228 and d = -33.9069 (taken from 
[94]). 

III. GROUND STATE PROPERTIES 

In Fig. 1 we present an overview of all the nuclei cov- 
ered in the present work. The chosen region ranges from 
the heaviest trans-actinides, well known from numerous 
experiments, up to beyond the neutron magic shell num- 
ber TV = 184 predicted by many theoretical models. The 
upper limit is provided by the vanishing of the two-proton 
separation energy. The latter approximates the proton 
drip line which is determined by the one-proton separa- 
tion energy with some correction due to the influence of 
the centrifugal barrier [95]. The region of neutron rich 
isotopes beyond the /3-stability line (indicated in Fig. 1 
by black squares) has been omitted as it is out of reach 
for the current experimental methods and can not be pro- 
duced in heavy ion collisions. As it can be seen in Fig. 1 
the experimentally known SHE's are located in the center 
of the investigated region. Since we consider the ground 
state properties as specially relevant for the understand- 
ing of the underlying physics we present these properties 
in the first subsection while the following subsection is 
entirely devoted to the study of the fission barriers. 



The theoretical approach discussed in the previous sec- 
tion has been applied to perform a systematic study of 
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FIG. 1: (Color online) Isotopes considered in this work. 



the properties of 160 super- heavy nuclei in the region 
150 < < 190, 100 < Z < 126. Earlier calculations 
in this region have been performed in the microscopic- 
macroscopic approach [32, 33, 39], in the HF plus BCS 
with the Skyrme force in the p-h channel and a monopole 
pairing force in the p-p channel [41, 42, 44-48] and in the 
full HFB approach by [49, 50, 96] with the Gogny force. 
The calculations by Berger et al. are rather similar to 
ours in the basic aspects. However, in our work the con- 
figuration space is larger, more appropriate for fission, 
and we allow more general HFB wave functions (with 
more simultaneously broken symmetries). Furthermore 
our study is rather detailed and systematic. 



A. Single particles energies 

To gain insight into the relevant physics of the SHE's 
we display in Fig. 2 the single particles energies (s.p.e.) 
for the nucleus ^™Hs in reflection symmetric conflgura- 
tions as a function of the quadrupole deformation ^2 and 
the quadrupole moment Q2- We have chosen this nucleus 
as a representative of the whole region for the following 
reasons. Since it is in the center of the area of the cal- 
culated nuclei its s.p.e. spectrum is characteristic for all 
she's, furthermore its proton and neutron pairing en- 
ergies are small and the extraction of the s.p.e. from 
the HFB calculations is more reliable. The s.p.e. are 
obtained as usual: After solving the HFB equations the 
one-body HF Hamiltonian is diagonalized in that basis. 
For the neutron system we find large spherical gaps at 
neutron numbers 164, 184 and 228; prolate gaps at 162 
and oblate gaps at 172 and 178. In nuclei with neutron 
number close to 160 the spherical minimum at 164 is over- 
whelmed by the prolate one at 162. For protons there are 
several gaps. The more relevant are the spherical ones for 
Z values 92, 114 and 126, the spherical-oblate one at 120 
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FIG. 2: (Color online) Single particle energies of the nucleus ^^''Hs as a function of the deformation parameter j32 (top x-axis) 
and the quadrupole moment Q2 (bottom x-axis). Continuous (dashed) lines represent positive (negative) parity states and the 
bullets the Fermi level. 



and the prolate ones at 104 and 108. As we will see be- 
low many nuclear properties can be understood just by 
looking at these s.p.e.'s. 



B. Deformations 

Let's now discuss some ground state properties start- 
ing with the deformation parameters /?„. Since most 
of the studies on these parameters have been done in 
the macroscopic-microscopic models we will use a def- 
inition of the in the spirit of these models. The 
deformation parameters /3„ entering in R{9) = _Ro[l + 
Sn=2 PnPn{cos6)\ are determined in such a way that the 
multipole moments, {Q\, A = 2 — 10), calculated with 
R{9) and the HFB self consistent ones do coincide. The 
deformation parameters of the calculated SHE's are 
given in the third column of Table I and they are visual- 
ized in Fig. 3a. The quadrupole deformation parameter 
/32 barely depends on the proton number Z and decreases 
slowly with increasing neutron number. For all Z values 
the nuclei with A'^ < 170 have a prolate deformation with 
132 > 0.25. In the neutron number region 170 < N < 182, 
the potential energy surfaces (PES) [51, 97] of these nu- 



clei present coexisting prolate and oblate minima. The 
quoted values for those nuclei are the ones corresponding 
to the deeper minimum. The prolate minimum is deeper 
in lighter nuclei with TV < 174 and Z < 118, while the 
oblate one is deeper in the heavier ones. With increasing 
neutron number the depth of the prolate well becomes 
smaller compared to the oblate one with the oblate min- 
imum becoming the lowest one around N=174 (see also 
Fig. 10). At the same time the absolute value of the 
quadrupole deformation parameter of the ground state 
decreases from j/^zl = 0.2 for A^ = 172 to |/32| = 0.05 for 
A^ = 182. The barrier between the two wells diminishes 
and finally disappears at A^ = 184, where the nuclei be- 
come spherical. We also observe in this region a weak 
Z dependency. For isotopes with A^ > 184 small pro- 
late deformation can be observed in the ground state. 
As compared with other calculations [98-101] of SHE's 
our /?2 are somewhat larger but this may be connected 
to the slightly different methods of computation of the 
deformation parameters /3„. In the presence of coexist- 
ing minima, one must be aware that when the restric- 
tion to axially symmetric shapes is released one of the 
two minima may not be a true minimum, but rather a 
saddle point [91]. A clear example is the case of the 
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FIG. 3: (Color online) Ground state deformation parameters P2, Ps, Pi, and Pe for super-heavy nuclei. Notice that the 
deformation scale is panel dependent. 



heavier Ds isotopes. In Fig. 10 one observes two min- 
ima (prolate and oblate) around ±5 b, but looking at 
triaxial calculations one finds that the prolate minimum 
is a saddle point. The gross behavior of the deformation 
parameter can be qualitatively understood by looking at 
Fig. 2. Here we observe that shell effects favor strong 
prolate deformations for the nuclei with 150 < N < 162 
and 100 < Z < 108, i.e., deformation driving (down- 
sloping) levels are being populated and spherical-driving 
(up-sloping) depopulated. In the region 164 < N < 170 
and 108 < Z < 116 neutron shell effects favor smaller 
prolate deformations while in the proton system small 
prolate as well as spherical shapes are favored. Around 
iV = 172 the oblate shapes are favored while larger N 
values prefer smaller oblaticity. This effect is reinforced 
by the soft oblate proton gap at Z = 120. The spheri- 
cal shape is the privileged one by neutron shell structure 
at N = 184. This coincides with the subsequent proton 
spherical gaps at Z = 114, Z = 120, and Z = 126 which 
cover most of the N = 184 isotones considered here. 

In columns 4 to 7 of Table I as well as in Fig. 3b,c,d 
we present the deformation parameters of higher multi- 
polarities. Most of the nuclei are reflection symmetric in 
the ground states and only for nuclei with N > 184 we 
obtain octupole deformed ground states independent of 



the proton number. MoUer et al. [98] also found non-zero 
deformations for the SHE's in this region. The odd 
multipolarities 13^, (ij and /Sg, not shown here, are differ- 
ent from zero only for those nuclei whose ground state 
is octupole deformed, however the numerical values are 
very small (e.g. jS^ values are less than 0.009). Concern- 
ing the even deformations higher than two, see Fig. 3c, d, 
they are different from zero only for nuclei with 7^ 0, 
i.e., for all nuclei with the exception of the few isotopes 
around N = 186. As we can observe there is a smaller 
dependence with Z than with N . Nuclei with N w 150 
have small negative hexadecapole deformations. For in- 
creasing iV, I/34I increases up to TV = 168 where it reaches 
the largest absolute value. Large negative hexadecapole 
moment together with large positive quadrupole mo- 
ment produces barrel-like shapes in the ground state 
of nuclei in this region. From here on |/34| decreases 
rather smoothly up to = 176 where it sharply de- 
creases to very small deformations and zero values at 
the largest neutron number analyzed. This sharp de- 
crease is associated to the prolate oblate shape transition 
that takes place at this neutron numbers. For A^-values 
150 < N < 176, i.e., the prolate part, there is a clear 
tendency with Z for a given N: with increasing proton 
number l/S^l gets larger. For the oblate-spherical part. 
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there is almost no dependence with Z . The behavior of 
the deformation parameter /3g is quite different from fi^. 
For 150 < N < 160 we obtain negative values decreasing 
in absolute value as A'' increases. Around A'' = 162, (3q 
becomes almost zero and for 164 < N < 174 we obtain 
increasing values of as increases. For N < 176 we 
obtain rather small /Je deformations. Our results for ^4 
and /Se are rather similar to those of MoUer et al [98] 
and Sobiczcwski et al [102]. Contrary to the (34 and (3^ 
parameters, /Jg is positive in the region 150 < N < 162. 
The nuclei with N = 164 and 166 have either zero or 



very small /?§ values. The nuclei with 168 < N < 17 A 
have mostly negative /?§ values. With the exception of 
the N = 178 isotopes that have positive deformations, 
all the nuclei with A'" > 176 have /Jg deformations al- 
most zero. We have to notice that the maximal absolute 
values of the deformation parameters decrease with in- 
creasing multipolarity indicating a decreasing relevance. 
This means that higher multipolarities can be omitted 
in the minimization of the energy of the ground state in 
non-self-consistent approximations. 



TABLE I: Selected properties of super-heavy nuclei: ground state defor- 
mation parameters, pairing energies, Qa, 2-micleon separation energies 
and half-lives for a emission and spontaneous fission calculated in HFB 
theory. The energies are given in MeV and half-lives in s. 
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FIG. 4: (Color online) Ground state proton and neutron pair- 
ing energies for super-heavy nuclei. 



C. Pairing energies 



of Fig. 4. The general behavior of these results can be 
easily understood by looking again at Fig. 2. We know 
that small pairing energies correspond to situations of low 
level density of the neutron s.p.e. In particular we expect 
zero neutron pairing energies at = 162, 178, 182 and 
184. Notice that the fact that for A — 162 the pairing 
energies are zero is consistent with the prolate deforma- 
tion of these nuclei and for A = 178 with the oblate one. 
For A = 182 and 184 we have spherical nuclei. Appro- 
priate neutron shell gaps can be found in Fig. 2 for the 
mentioned cases. For the same reason large pairing ener- 
gies are associated with high level density, for example, 
we find large pairing energies for A = 158, 168 and 190. 
We observe that the general pattern, qualitatively, does 
not depend on the proton number. The proton pairing 
energies are given in the ninth column of Table I and the 
lower panel of Fig. 4. In this case the general pattern 
looks more complicated. For the same proton number 
some nuclei are prolate, some oblate and some spherical 
and therefore only part of them fit to the proton shell 
gaps in Fig. 2. Again the small pairing energies found at 
Z values 104, 108, 116, 120 and 126 have to do with the 
low proton level density found in Fig. 2 for those pro- 
ton number. Concerning the behavior for the different 
isotopes, for example for Z = 112, from A = 158 up 
to A = 174 we have proton pairing energies around —9 
MeV and from A = 176 on around —20 MeV. The rea- 
son for this change is that in the first interval the nuclei 
are prolate and in the second either oblate or spherical. 
Similar arguments apply for other cases. It is interesting 
to notice that with the exception of the few nuclei with 
A > 188 none of the nuclei studied has absolute values 
of the neutron pairing energies larger than 9 MeV. On 
the other hand, many of the analyzed nuclei have proton 
pairing energies much larger than 10 MeV. 



In the HFB approach [81], the pairing energy is given by 
= -irr(AK*), (9) 

with 

Afejfe^ = 2 X! ^kik2k3kiHk3ki, (10) 

the pairing field and 

Kk,M = {HFB\ck,CkAHFB), (11) 

the pairing tensor. 

The neutron pairing energies of each nucleus are given 
in the eighth column of Table I and in the upper panel 



D. Qa values 



In Fig. 5 the Qa values are plotted as a function of the 
neutron number [see Eq. (8)]. The exact numerical val- 
ues are also given in column 10 of Table I. The patterns 
displayed by the different isotopes are easily understood 
just looking at Eq. (8): minima appear when the mother 
nucleus are more bound than the average and maxima 
correspond to more bound daughter nuclei. Thus we ob- 
serve minima corresponding to the neutron numbers for 
which energy gaps appear in the single particle energies 
of Fig. 2, namely 162, 178 and 182-184. For the A values 
164 and 172 there is no structure because the energy gaps 
that one finds for these neutron numbers are not large 
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FIG. 5: (Color online) Qa values, theoretical values marked 
by dots are compared with experimental data [10, 93]. The 
different lines correspond to the indicated isotope 



enough to provide energy minima at those deformations, 
see column three of Table I. Experimental values [10, 93] 
for some isotopes of the nuclei Fm, No, Rf, Sg, Hs, Ds, 
114 116 and 120 are also displayed in the Figure. The 
agreement between theory and experiment is very sat- 
isfactory, in most of the cases we obtain a quantitative 
agreement and for the others at least the tendency is the 
right one. Concerning the proton dependence we observe 
a similar situation, we find large energy spacing for the 
lines whose Z numbers correspond to energy gaps in the 
single particle diagram of Fig. 2, for example, Hs, 116, 
120, etc. 



E. Two-nucleon separation energies 
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FIG. 6: (Color online) Top panel: Two neutron separation 
energies as a function of the neutron number. Lower panel: 
Two proton separation energies as a function of the proton 
number. The lines represent isotonic chains, the one at the 
lower left corner corresponds to N=150, and the one closer 
to the upper right corner to A=190. Experimental data are 
taken from Ref. [93]. 




In Fig. 6, upper panel, we present the 2-neutron sep- 
aration energy as a function of the neutron number (see 
also column 11 of Table I). As expected we obtain a de- 
creasing behavior of S2N with increasing neutron number 
since we get closer to the neutron drip line. The smooth 
decline of S2N is only disturbed at the neutron numbers 
corresponding to the single particle shell gaps 162, 184, 
etc. The two-proton separation energies are shown, in 
the lower panel of Fig. 6, as a function of the proton 
number (see also column 12 of Table I). The general be- 
havior of decreasing 6*2 p for a given isotonic chain with 
growing Z illustrates the fact that we are getting closer to 
the proton drip line. The fact that the S2P presents less 
structure than the S2N is obviously related with the fact 
that, in the region of interest, the shell gaps in Fig. 2 
are smaller for protons than for neutrons. Fig. 6 also 
includes the available experimental values for some Fm 
and No isotopes. We observe that in the case of S'2Ar the 



theoretical values are a around 1 MeV smaller as com- 
pared to experiment while in the case of S2P they are 
around 1 MeV larger. In both cases, however, the trend 
is correctly described. This is the well known binding 
energy drift that takes place with the DIS parametriza- 
tion of the Gogny force for most isotopic chains. To cor- 
rect this drift Chappert et al. [103] have recently done a 
new parametrization of the Gogny force, the DIN, which 
reduces considerably the drift but otherwise keeps the 
quality of the DIS parametrization. In Table II we have 
included calculations with the DIN parametrization and 
we observe a considerable improvement in the agreement. 
We also include in the table the results of calculations 
with the DIM parametrization of the Gogny force ob- 
tained by including beyond mean field effects in the fit 
[104], which also do not present the mentioned drift. 

The DIN and DIM parametrizations also improve 
slightly the agreement of Qa with the experiment (see 
Table II), but this effect is much less pronounced than 
in the case of the separation energies as the influence of 
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TABLE II: Twonucleon separation energies and Qa in Fm 
isotopes calculated with various Gogny forces compared with 
experimental data [93] (units are MeV). 





N 


DIS 


DIN 


DIM 


EXP 




150 


13.20 


13.86 


14.00 


13.98 




152 


12.32 


13.01 


13.14 


13.40 




154 


11.35 


12.13 


12.23 


12.06 




156 


10.77 


11.64 


11.73 


11.56 




150 


8.86 


7.80 


7.62 


7.74 




152 


9.93 


8.79 


8.59 


8.93 


S2P 


154 


10.87 


9.67 


9.47 


9.71 




156 


11.74 


10.50 


10.31 


10.43 




150 


7.35 


7.66 


7.70 


7.56 




152 


7.12 


7.49 


7.54 


7.15 


Qa. 


154 


7.02 


7.38 


7.48 


7.31 




156 


6.66 


6.99 


7.10 


7.03 



protons and neutrons cancel each other. 



IV. SPONTANEOUS FISSION 
A. Symmetry breaking effects in the fission paths. 

The large size basis used in these calculations and the 
large number of nuclei studied prevents to make a system- 
atic study with totally unrestricted symmetry-breaking 
wave functions. To clarify the effects of these restrictions 
we show in Fig. 7, as an example, the behavior of several 
relevant quantities along the fission path for the nuclei 
274Hs (left panel), ^^Sds (middle panel) and ^^'^Cn (right 
panel) in different approximations. The panels (d), in 
the middle row of Fig. 7, display three different fission 
paths for each nucleus corresponding to the following 
constraints on the wave functions: 1.- The Axially Sym- 
metric (7 = 0) and Reflexion Symmetric {Q3 = 0) fission 
path (continuous line), which we shall call in the follow- 
ing AS-RS path. 2.- The Axially Symmetric and Non- 
Reflexion Symmetric {Q3 ^ 0) path (long-dashed lines), 
in the following AS-NRS path. 3.- The Non- Axially Sym- 
metric (7 7^ 0) and Reflexion Symmetric path (short- 
dashed lines), in the following NAS-RS path. 

In the ground state these nuclei are well quadrupole 
deformed, Q2 ~ 15 b, see panels (d), and in the way 
to fission they have to tunnel through a barrier of sev- 
eral MeV height. The barriers represent the potential 
energy needed to deform the nucleus. It is related, there- 
fore, with the single particle levels around the Fermi sur- 
face available at the corresponding deformation along the 
path. For the three nuclei the paths in the AS-RS case 
present two hump barriers where the height of the sec- 
ond barrier decreases with the mass number while the 
first one remains more or less constant for ^^^Hs, ^^*Ds 
and increases for ^*^Cn. The width of the barrier in this 
approach is also similar for the three nuclei. The origin 
of the two humps can be easily understood looking at the 



single particle energies of Fig. 2. If we follow the neutron 
Fermi level to the prolate side we find two regions with 
a clearly dcivclopcd low level density. The first one at 
(52 ~ 15 b with the shell gaps N=162 and Z=108 corre- 
sponds to the ground state minimum. On the way from 
this point to larger deformations we find a larger level 
density region which corresponds to the energy increase 
of the first barrier. Behind that we arrive at the second 
region at Q2 ~ 32 b, which corresponds to the super de- 
formed minimum. At even larger deformations one finds 
a high level density region that provides the energy rise 
of the second barrier. Finally above 40 b, we observe 
intruder states of a high lying I'i^- orbital which could 
be an indication of the scission point. Incidentally, since 
the neutron shell gap is at N=162 and the proton one at 
Z=112 the super-deformed minimum is deeper in ^®^Cn 
than in the lighter nuclei displayed in Fig. 7. 

This behavior changes remarkably in the AS-NRS ap- 
proach. One can follow in panels (f) the portion of the 
paths where a lower solution with Q3 7^ is found. As 
one can see in panels (f) the octupole degree of free- 
dom plays an important role for quadrupole deformations 
starting around Q2 ~ 30 — 35 b, i.e., close to the super 
deformed minima. As a matter of fact in this channel 
the second hump of the barriers diminishes strongly and 
the paths in this region look like the continuation of the 
first hump. The AS-NRS path reaches up to large defor- 
mations, Q2 « 108 b for 274hs, « 77 b for ^^SDs and 
decreases to Q2 « 67 b for ^^^Cn. 

We can understand the onset of octupole deforma- 
tion looking again at Fig. 2. Here we observe that the 
K = 1/2 and K = 3/2 levels of the i^kjj; shell cross 
the Fermi surface around 26 b, while the K = 5/2 level 
crosses around 35 b. Interestingly, at these Q2 values 
the K = 5/2 level of the i^h.^- orbit also approaches the 
Fermi surface. In Fig. 2 and at zero deformations the 

vkiz subshell lies at about 0.4 MeV and the vh^ at 

2 ^ 

about —10.4 MeV. The ki? and subshells interact 

2 

strongly through the AL = A J = 3, i.e., the octupole 
interaction. That means if we allow for reflexion symme- 
try breaking we can increase the quadrupole deformation 
at a lower energy cost [105]. We can observe in Fig. 7 
that around these values the AS-NRS fission paths get 
lower in energy than the AS-RS ones. In the (g) panels 
we can follow the behavior of the hexadecapole moment 
along the fission path. In the ground state it is close to 
zero in both approximations and from this point on it 
grows linearly with Q2 in the AS-NRS approach. In the 
AS-RS approach, however, first it increases linearly up 
to the scission point where a kink is observed followed by 
a linear increase. 

In the NAS-RS path we can observe the effect of the tri- 
axial shapes along the fission paths. In panels (e) we can 
see the two portions of the trajectory where triaxial solu- 
tions are found. The first one, close to the ground state, 
spans a smaller interval of Q2 values than the second one 
and does not have a large impact on the energy. The sec- 
ond one, as in the AS-NRS case, is relevant around the 
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Q2 [b] Q2 [b] Q2 [b] 

FIG. 7: (Color online) Behavior of several magnitudes for the nuclei ^^''Hs, ^^*Ds and ^^^Cn along different fission paths. 
Neutron and proton pairing energies Ep (a) and Ep (b) are in MeV, (c) inertia parameter B in lQ~^h~^ MeV~^ , (d) HFB 
potential energy in MeV, (e) triaxial deformation parameter 7 in degrees, (f) octupole moment Qz in b^^^, (g) hexadecapole 
moment Q4 in , (h) sub-integral function of action integral dS/dQ2 in . 



second hump causing a significant lowering of its height, 
i.e., in this case we still have to deal with two humped 
barriers. The width of the barrier, at variance with the 
AS-NRS case, is more or less like the AS-RS one. 

Though the shape of the barriers is very relevant to cal- 
culate lifetimes one has to consider, however, that other 
quantities entering in the corresponding formula, Eqs. (1) 
and (2), do play an important role. A relevant parameter 
is the collective quadrupole inertia B{Q2), since mass pa- 
rameters are strongly influenced by pairing correlations 
and these itself by the single particle level density which 
we expect to vary along the different fission paths. In 
panels (a) and (b) of Fig. 7 we display the neutron and 
protons pairing energies, respectively. We indeed observe 
big differences in both of them along the portions of the 
paths where the symmetry breaking takes place. In par- 
ticular we observe that the AS-RS solutions provide al- 
ways the largest pairing energies. As we can observe in 
panels (c), where the i?(Q2)'s are plotted in the differ- 
ent approaches, in the relevant parts the AS-RS provides 
the smaller masses followed by the AS-NRS ones. This 
implies, that not necessarily the smallest fission barri- 
ers provide the shortest lifetimes [29]. Since the action 



S in Eq. (2) can be seen as a line integral of the func- 
tion dS{Q2)ldQ2 = {2B{Q2){V{Q2) - Eofl^, the area 
of the surface delimited by this line and the X-axis pro- 
vides the value of the action 5*. In panel (h) we display 
dS{Q2)/dQ2 in the three approaches. For the nucleus 
^^''Hs wc find that though the NAS-RS fission path has 
a smaller fission barrier than the AS-RS, the actions S 
for both approaches are very close. The actual values 
of S are 26.69 (AS-RS) and 26.49 (NAS-RS), while in 
the AS-NRS case we obtain a much larger value, namely 
41.66. In the nucleus ^^*Ds, the fission paths alone would 
predict that the AS-NRS and the NAS-RS approaches 
would provide much shorter lifetimes than the AS-RS 
one. However, in panel (h) one finds that the three areas 
look rather similar. Actually the precise numbers 25.88 
(AS-RS), 26.97 (AS-NRS) and 26.49 (NAS-RS) show this 
to be the case. Lastly, for ^^-^Cn, the prediction of the 
fission paths is more or less in accordance with the one of 
panel (h) and the actual numbers 23.34 (AS-RS), 18.32 
(AS-NRS) and 22.52 (NAS-RS) corroborate that. We 
can conclude that the restriction to axially symmetric 
paths is, in general, a good approximation, though as we 
we will see later one can find some exceptions. 
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The low and short NRS barrier in nuclei with A'^ > 170 
makes the most probable fission through octupole de- 
formed shapes. In these nuclei we expect to find asym- 
metric mass distribution of fission fragments [106]. 

In the one dimensional fission paths plotted in Fig. 7 
wc find crossings between the two paths giving the im- 
pression that one could switch from one path to the other 
without further problem. However, if we look in a higher 
dimensional plot one can see that this is not the case. 
To illustrate this point we have drawn in Fig. 8 potential 
energy contour lines versus the quadrupole, Q2, and the 
octupole moment, Q3, for the nuclei ^^'*Hs and ^^^Cn. 
In this figure we can follow the AS-RS and the AS-NRS 
paths of Fig. 7 for the respective nucleus. The AS-RS 
path corresponds to Q3 = and goes along the X-axis 
and the AS-NRS one goes along the bullets. It is interest- 
ing to see how the self-consistent path goes along a valley 
in both nuclei. We can also see that no alternative paths 
are present. In the ^''^Hs case we find that, at Q2 = 50 
b where both paths seem to cross in Fig. 7(d), in reality 
both paths are separated by a 4 to 5 MeV high barrier. 



B. Fission barriers in the axially symmetric 
approaches. 

In the following to perform a systematic description of 
the fission barriers of the 160 SHE's, we restrict ourselves 
to the axial approximation in which we have performed 
two kinds of calculations, namely the reflexion symmet- 
ric, AS-RS, and the non-reflexion symmetric, AS-NRS. 
All the fission barriers are presented in the series of Fig- 
ures 9-13. 

In Fig. 9 we present the fission barriers for the isotopes 
of the elements Fm, No and Rf for quadrupole values 
from —20 b up to 80 b (continuous line for AS-RS and 
dashed for AS-NRS). Wc first discuss the AS-RS rcsuhs. 
In panel (a) we present the Fm results for neutron num- 
ber 150 up to 166. All isotopes present a well prolate de- 
formed minimum around 15 b. In addition in the lighter 
isotopes a shallow superdeformed (SD) minimum appears 
around 50 b, at iV = 156 we find a very fiat minimum and 
for the heavier isotopes no SD minimum is found. The 
common characteristic of these nuclei is the presence of 
a big broad barrier. For N = 150 the barrier is centered 
at Q2 = 30 b, has a height of about 12 MeV and a width 
of 18 b. With increasing neutron number the center of 
the barrier shifts to larger deformations and the height 
diminishes. For N — 166 the center is around Q2 — 38 h 
and the height is about 8 MeV. In the heavier isotopes we 
find some structure in the first barrier, namely the devel- 
opment of a shoulder around Q2 = 27 b with increasing 
neutron number. The presence of a SD minimum in the 
lighter isotopes drives the existence of a second barrier. 
Since the minima are rather shallow the second barriers 
are broad but not high. These properties will contribute 
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FIG. 8: (Color online) The PES of ^^Hs as a function of 
quadrupole and octupole moments. The energy origin has 
been set at the energy minimum. 

in general to a tendency of shorter lifetimes with increas- 
ing neutron number, though the particular behavior must 
be analyzed case by case. 

The fission paths for the No isotopes are shown in panel 
(b) for the same neutron numbers as the Fm case. The 
structure of the first barriers look roughly like the ones in 
Fm. An important point is the fact that the SD minimum 
is somewhat lower in energy than in the Fm case and its 
role is therefore much less relevant. The results for the 
Rf isotopes are displayed in panel (c). The tendency ob- 
served in the No isotopes is reinforced, the SD minimum 
gets even deeper and the second barrier disappears for 
most isotopes. This fact will provide in general shorter 
lifetimes for Rf than for No and for No shorter than for 
Fm. These are general tendencies but since the lifetimes 
are very sensitive to small energy differences along the fis- 
sion path, to make quantitative predictions calculations 
involving also collective inertia have to be performed. 

Let us now describe the AS-NRS results. For all 
isotopes of the three elements, the fission barriers are 
much larger along the non-refiexion symmetric paths 
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FIG. 9: (Color online) Fission barriers for the nuclei Fm, No and Rf along the AS-RS (continuous lines) and the AS-NRS path 
(dashed lines). The ground state energy of each isotope has been set to zero and the barriers of the different isotopes have 
been shifted in 10 MeV in the axis of coordinates for clarity reasons. 



than along the reflexion symmetric path. The lifetimes 
are therefore considerably longer in the AS-NRS path 
than in the AS-RS one. The fission in the AS-NRS mode 
is completely impossible. 

An exhaustive discussion of the fission barriers of the 
nuclei 254Fni, ^^^Fm, ^^^Fm, ^^^No, and ^eoRf can be 
found in Ref . [54] . In these nuclei the "elongated fission" 
mode can be observed. This mode is connected to oc- 
tupole deformed fission paths which start at Q2 > 50 b. 
We have not observed the "elongated fission" in heavier 
nuclei therefore we will not discuss it here. 

In Fig. 10 we present the paths for the nuclei of Sg, 
Hs and Ds elements. In panel (a) we display the results 
from — 150 up to A = 174 for the Sg isotopes. These 
nuclei have a trend somewhat different from the preced- 
ing ones: the SD minimum does not play a relevant role 
in the fission process since it is always deeper than the 
ground state and as a matter of fact if the SD minimum 
were not that shallow it would be the ground state for the 
lightest Sg isotopes. Furthermore the barriers get flatter, 
most of them being lower than 10 MeV. Apart from this 



feature the situation for the isotopes A^ = 150 — 166 is 
similar to the nuclei Fm, No and Rf, i.e., the ground 
states are prolate deformed (Q2 = 15 b) with large nega- 
tive deformation energies and high barriers at larger Q2 
values. For A^ — 168—174 the situation changes very fast, 
and the following properties get reinforced as the neutron 
number increases: the ground states get less deformed, 
the shoulders around Q2 = 30 b get deeper and become 
real minima and as a consequence the original barrier 
becomes a two humped one. The height of the spherical 
maximum decreases and a coexisting oblate minimum 
develops at Q2 — —10 b. In particular for N = 172 the 
oblate and prolate minima are degenerated. We there- 
fore expect a strong reduction in the fission lifetime of 
heavier isotopes as compared to the lighter ones. 

In panel (b) the A^ — 150 up to A^ = 186 Hs isotopes 
are shown. In the lighter isotopes the fission barriers 
are flatter and the SD minima deeper than in the corre- 
sponding Sg isotones. The effect of the SD minima on 
the flssion process is much smaller. Furthermore the sec- 
ond flssion barrier develops much earlier. Since the single 



16 




-20 20 40 60 -20 20 40 60 -20 20 40 60 



Q2 [b] Q2 [b] Q2 [b] 

FIG. 10: (Color online) Fission barriers for the Sg, Hs and Ds isotopes along the AS-RS (continuous lines) and the AS-NRS 
path (dashed lines). The green short-dashed lines around the ground state for A > 184 correspond to NRS solutions. The 
ground state energy of each isotope has been set to zero and the barriers of the different isotopes have been shifted in 10 MeV 
in the axis of coordinates for clarity reasons. 
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FIG. 11: (Color online) Same as Fig. 10, but for the Cn, Z = 114 and 116 isotopes. 



particle energies of Fig. 2 are for ^^"^Hs we can confront 
the general discussion of Sect. IV A with the shape of the 
fission path in this nucleus. Aside from these facts the Hs 
isotopes from = 150 up to = 174 behave to a large 
extend like the corresponding isotones of the Sg isotopes. 
For ./V = 176 the oblate minimum becomes the ground 
state and the two humps of the fission barrier are very 



similar in size and height (about 5 MeV) . With increasing 
number of neutrons, the prolate minimum shifts towards 
Q2 = and the spherical minimum becomes the ground 
state at the A^ — 184 shell closure. As a consequence 
the two humps separate from each other, the outer one 
shifting to larger Q2 values and the inner one to smaller 
ones. At the same time the inner barrier gets bigger and 
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the outer smaller. From this behavior one expects the 
lifetimes for fission to be smaller around neutron number 
170-174. 

The nucleus ^^''Hs is octupole deformed in its ground 
state, see Fig. 3. In the self-consistent (52-constrained 
calculations this nucleus remains weakly octupole de- 
formed up to (^2 = 10 b where it turns reflexion sym- 
metric (see short dashed part of the fission path). Since 
the fission fragments at the scission point are character- 
istic of a refiexion symmetric fission we still denote this 
mode as RS fission. The same situation is found with the 

> 184 isotoncs for the heavier SHE. The paths along 
Qa close to the ground state are plotted with short 
dash lines in the corresponding figures. 

The Ds results for the isotopes N = 154 up to TV = 188 
are displayed in panel (c). The main characteristics of 
these fission paths are the following: For the lighter iso- 
topes the fission barriers are flatter than for the cor- 
responding isotones in Hs and Sg. For the light and 
medium mass isotopes, for a given isotone number, we 
find an increase of the first barrier moving from Sg to 
Hs and from this to Ds. The opposite effect is observed 
for the second barrier, in particular for the very heavy 
isotopes this barrier disappears at the highest neutron 
number studied. 

Concerning the AS-NRS results of Fig. 10 and for the 
Sg isotopes, we find that though the AS-NRS fission bar- 
riers are smaller than the AS-RS ones for medium Q2 
values, for larger ones they are much higher and as a 
result the non-reflexion symmetric fission paths are not 
favored as compared with the reflexion symmetric ones. 
For the Hs isotopes the same behavior as the Sg isotopes 
is observed: for lighter isotopes up-bending tails of the 
fission paths make the AS-NRS paths very imfavored, 
but for medium-heavy N values the tails bend down and 
though the path is longer the second barrier is smaller for 
the AS-NRS path than for the AS-RS one, i.e., around 
A'' = 170 and above the AS-NRS path may compete with 
the AS-RS one (cf. Fig. 8). For the heaviest Hs isotopes 
in the AS-RS approach the second barrier decreases con- 
siderably and the AS-NRS barrier becomes longer than 
the other one. For the Ds isotopes the AS-NRS path 
is even more favorable because the down-bending ten- 
dency gets reinforced and we have some AS-NRS paths 
which are clearly favored, for example for the nuclei with 
iV = 170 — 176. For larger N values the vanishing of the 
second barrier in the AS-RS case again favors this ap- 
proach. Notice that in general octupole effects set in for 
larger Q2 values as the neutron number increases. 

In Fig. 11 we present the fission paths for the isotopes 
of the elements Cn and Z = 114 and 116. As before we 
first discuss the AS-RS paths. In panel (a), where the 
isotopes A'' = 160 — 188 of Cn are shown, we observe that 
compared with the corresponding isotones discussed be- 
fore, the second hump of the barriers is lower and that 
the slopes of the tails of the fission paths are more pro- 
nounced. These facts point to shorter lifetimes of the Cn 
isotopes as compared with lighter isotones. In panel (b) 



of the figure the Z = llA isotopes are shown. The main 
difference with respect to the Cn isotopes is the increase 
of the first hump and the decrease of the second one. In 
the isotopes with one hump barriers these are higher and 
a bit broader than for the corresponding isotones in Cn. 
Altogether, it seems that, in general, the lifetimes of the 
Z = 114 isotopes will be somewhat longer than the one 
at the corresponding Z = 112 isotones. In panel (c) we 
display the Z = 116 isotopes. Here the same trend as 
in the previous nuclei is observed: a reinforcement of the 
tendency to increase the first hump of the double humped 
barriers and in the case of only one hump an increase of 
this. 

Concerning the AS-NRS results for Cn and the Z = 
114 and 116 elements, we find that the onset of octupolar- 
ity is energetically favored after the level crossing of the 
"higher shells" and with the exception of the lightest iso- 
topes the barriers are much smaller in the AS-NRS path. 
This is due to the disappearance of the second hump of 
the barrier, i.e., in the AS-NRS path we have only one- 
humped barriers. We expect therefore a shortening of 
the fission lifetimes along these paths. 

In Fig. 12 the results for the Z = 118, 120 and 122 iso- 
topes is shown. The same tendency as before is observed 
in the AS-RS calculations, larger first barriers as Z in- 
creases. The role of the octupole degree of freedom is also 
relevant and all AS-NRS fission paths do have smaller fis- 
sion barriers. The same comments do also apply to the 
Z = 124 and 126 isotopes in Fig. 13. 



V. HALF-LIVES OF SHE'S 

One of the key issues in the theoretical description 

is the prediction of the decay modes and the half-lives 
of the she's. The agreement between the experimen- 
tal data and the theoretical predictions states not only 
the quality of both of them but it may be also a cri- 
terion for the identification of a particular isotope syn- 
thesized in the corresponding fusion reaction. Moreover, 
since the contemporary experimental techniques do not 
allow the detection of nuclei with half-lives shorter than 
T = 10 /is one must estimate which isotopes can survive 
long enough to be detected and the way in which they 
disintegrates. Therefore, the half-lives of the two main 
competing processes, spontaneous fission and a emission, 
should be evaluated. The shortest half-life determines 
the dominant decay channel and the total half-life. If the 
branching ratio between two modes is equal to 50 % the 
logarithm of the total half-life would be smaller than the 
shorter partial half-life by not more than 0.3. Differences 
in half-lives between two modes of one order of magni- 
tude lead to a logarithm of the total half-life only 0.05 
smaller than the logarithm of the half-life of the fastest 
decay. 

The half-lives for a decay and spontaneous fission cal- 
culated in the HFB theory are collected in the last two 
columns of Table I. For an easier analysis these data are 
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FIG. 12: (Color online) Same as Fig. 10, but for the Z = 118, 120 and 122 isotopes. 



also presented in Fig. 14 where the isotopic chain of each 
element is shown in a separate panel. In some nuclei two 
fission half-lives can be calculated along paths leading to 
distinct fragment mass asymmetry. Both solutions are 
depicted in Fig. 14 and discussed below, however, in Ta- 
ble I only the shorter half-life of the dominant mode is 
given. The available experimental data are also plotted 
in Fig. 14. 

A first look into the panels of Fig. 14 reveals similar 
tendencies for the different isotopic chains. It is easy to 
distinguish the intervals of neutron number where com- 
mon features are characteristic for many elements despite 
of differences in the absolute values. Therefore we will 
discuss our results collected in groups of similar neutron 
number starting from the lighter ones. 

As we have seen in the discussion of the precedent sec- 
tion NRS effects influence the fission paths, and thereby 
the fission half-lives, in three Q2 regions: 1.- For small 
Q2 values, Q2 < 20 b, they affect nuclei with octupole 
deformed ground states, i.e. nuclei with N > 184. 2.- 
Starting at medium {Q2 > 20 b ) and extending up to 
large Q2 values, these effects are present in all nuclei. 3.- 
For larger Q2 values, Q2 > 50 b, they appear in the light 



isotopes of the elements Fm, No, Rf and Sg. According 
to these effects we are using three symbols in Fig. 14 for 
the fission half-lives depending on the paths used in the 
calculations : AS-RS (triangle-up) includes the genuine 
AS-RS ones plus those with octupole effects close to the 
ground state. The reason to include the latter under a 
"AS-RS" denomination is that at ~ 20 b the nucleus 
takes AS-RS shape and the fission takes place exactly in 
the same way as in the pure AS-RS case with symmet- 
ric fragment mass distribution. AS-NRS (triangle-down) 
denotes the "real" NRS paths leading to fission with frag- 
ments of different masses, which correspond to the long 
dashed lines in the fission barrier plots. AS-RS/NRS 
(star) labels the paths where the first part is AS-RS and 
only for > 50 b one follows the NAS branch. 



A. The region 150 < AT < 162 

The first characteristic region covers the lighter iso- 
topes with 150 < N < 162, clearly delimited by the 
N = 162 deformed "shell closure" . In this range of N 
the half-lives for a emission increase monotonically with 
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the neutron number in all elements from Fm to Z = 114. 
It can be also noticed that an increase of the proton num- 
ber by two units leads to the decrease of the a half-lives 
by two or three orders of magnitude. As we will also 
see for the other regions, this tendency is a direct con- 
sequence of the calculated values of Qa shown in Fig. 5 
and is consistent with well-known properties of a decay 
in heavy nuclei. 

The three lightest Fm and No isotopes are a specific 
group of nuclei in which the HFB calculations predict rel- 
atively long fission half-lives. The small fission probabil- 
ity is the consequence of the shape of the fission barrier in 
these nuclei which extends to large quadrupole deforma- 
tion. In the AS-RS path the second barrier extends up to 
120 b providing very long half-lives. A somewhat shorter 
half-life is obtained in the AS-NRS approach where the 
barrier is constructed from two humps: the first a RS and 
the second a NRS one starting at Q2 ~ 60 b (Fig. 9). In 
the slightly heavier isotopes the second NRS barrier also 
can be noticed in the PES (Fig. 10), but in contrast to 
the group of lighter Fm and No nuclei the second, su- 
per deformed minimum has an energy below the ground 
state. Thus the second hump does not affect the barrier 
tunneling but it governs the fragment mass asymmetry. 
In these nuclei the reduction of the width of the barrier 
leads to a decrease of the half-lives by a few orders of 
magnitude what can be seen in the Fm and No isotopes 
with N > 156 as well as the Rf and the Sg isotopes. All 
nuclei with the second NRS hump leading to the so-called 
"elongated fission" [54] (i.e. Fm and No with N < 158 as 



well as Rf and Sg with N < 156) are marked in Fig. 14 
by the blue stars. 

In the heavier isotopes of the 150 < N < 162 interval, 
where only the single RS barrier remains, one observes 
a rise of the fission half-lives with the neutron number. 
The slope of this trend changes from an almost flat de- 
pendence in Fm to a very steep one in Sg, Ds, Cn, and 
Z = 114. This trend is caused by a broadening of the bar- 
rier that eventually becomes a two humped one (Figs. 9, 
10, 11). 

The local maximum found in all elements in the partial 
half-lives at iV = 162 with respect to a and fission decays 
indicates the special character of this neutron number in 
the chart of nuclides, see also Fig. 2. These isotones are 
the most stable in the close vicinity, hence, in many pa- 
pers N = 162 is called "deformed magic number" . This 
name stresses the significant difference from the classical 
magic nuclei which are spherical in their ground states. 

Many experimental data coming from the "cold fu- 
sion" reactions are available in the region of isotopes 
with N < 162 [9, 107]. The agreement of the a decay 
half- lives with the theoretical predictions is noticeable al- 
though some discrepancies are observed. Most calculated 
fission half-lives overestimate the experimental data. The 
agreement would be better if triaxial effects were taken 
into account in the saddle point of the first barrier. The 
consideration of the 7-deformation in the calculations of 
nuclei in the Fm region reduces the barrier heights by 
around 2 MeV decreasing thereby the theoretical fission 
half- lives by around 2 orders of magnitude [54] . At vari- 
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FIG. 14: (Color online) The spontaneous fission and the a emission half-lives in SHE's. Experimental data are taken from 
Refs. [9, 10, 107] 



ance with the discussion of Sect. IV A in these nuclei this 
reduction is not fully compensated by an increase of the 
inertia parameters. Therefore discrepancies between the 
theory and the experiment diminishes considerably. It 
should be pointed out that the fastest decay is properly 
predicted in all cases. Furthermore, when the experimen- 



tal data for a emission and fission provide comparable de- 
cay probability for both processes the theoretical predic- 
tions are also similar for both partial half-lives. Almost 
all nuclei in this region do have half-lives long enough to 
be considered as of experimental interest. 
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FIG. 15; (Color online) The dominant decay modes ol SHE's. 
Logarithm of the shortest half-lives (in s) is also indicated. 



B. The region 164 < iV < 178 

The next region covers nuclei with neutron number 
from iV = 164 to iV = 178. At = 164 a kink of 2- 
4 orders of magnitude in the a decay half-lives can be 
observed corresponding to a local maximum in Qq,, see 
Fig. 5. The increase of the neutron number for a given 
isotope causes a linear growth in the a half-life following 
the tendency already observed in the lighter nuclei. The 
half-life rise up to the neutron number TV = 178 where the 
locally longest living isotopes are found, see also Fig. 5. 
They reach values from T = 0.1 s in Z = 118 to T = 
10^^ s in Hs. Again, the half-lives calculated for heavier 
isotones are smaller and, consequently, the a emission 
becomes the most probable decay channel in the proton 
rich nuclei. 

The AS-RS fission half-lives for nuclei with 164 < 
N < 178 behave completely different than those for 
150 < N < 162. After passing the maximum at A^ = 162 
they decrease up to the local minimum at A^ = 170. The 
drop is very steep for the lighter elements (Rf, Sg, Hs) 
and gentler for the heavier ones (Ds, Cn), whereas for the 
isotopes oi Z = 114 almost no change is observed. This 
behavior can be explained by the decreasing potential 
energy along almost the whole barrier, associated with a 
diminishing of the saddle point energy and the clear de- 
velopment of the second minimum (Figs. 9, 10, 11). The 
inverse trend of increasing the potential energy at small 
deformations is noticed in heavier elements, starting from 
Cn (Fig. 11). The first RS hump of the barrier grows up 
with increasing mass of a nucleus. In the heavy isotopes 
with A > 280 it is higher than the second RS barrier. 

The half-lives calculated for 162 < A^ < 170 along the 
NRS path diminish more rapidly with the increase on A^ 
than in the RS mode. This is induced by the changes that 
take place in the NRS barrier, namely the fast decrease 
of its height and, more important, the narrowing of its 
width, see Figs. 10, 11. In nuclei with N > 170 the NRS 



fission barrier allows to avoid the second hump of the RS 
barrier making the asymmetric fission the most probable 
mode. Very short fission half-lives, even below T = 1 ms, 
can be found in the Hs, Ds, Cn and Z = 114 isotopes, 
where the NRS fission is the dominant decay mode being 
even faster than the a emission. In Z — 114 the NRS 
fission has half-lives comparable to the a decay. 

The available experimental data around A^ = 172 ob- 
tained from the "hot fusion" experiments [8, 9], fit per- 
fectly to the theoretical predictions. In Cn two spon- 
taneous fission half-lives correspond to the prediction of 
the NRS fission mode. In Z = 114, A^ = 172 the ob- 
served 50% branching ratio is very well reproduced. In 
the isotope with two more neutrons the detected a de- 
cay is predicted with only slightly longer half-life than in 
the dominant fission channel. Finally, in Z = 116 and 
Z = lis the dominant a radioactivity is properly pre- 
dicted by the theoretical analysis with a good estimation 
of the half-lives. 

While approaching A^ = 178 the fission half- lives in- 
crease with a slope which grows with the proton number. 
This effect is governed by two trends observed in the 
evolution of the PES's. The first one is the aforemen- 
tioned growth of the first hump of the barrier which is 
the highest one in almost all nuclei with A^ > 170. The 
other factor is the lowering of the energy of the oblate 
minimum which becomes the ground state in the nuclei 
around A^ = 178 (cf. discussion of ground-state defor- 
mations in Sec. HI). The shift of the ground state from 
the prolate to the oblate minimum, see also Fig. 2, gives 
an additional increase of the barrier height up to over 1 
MeV. The energy difference between prolate and oblate 
minima has its largest value in the A^ = 178 isotones 
what, added to the large first barrier, produces extremely 
long half-lives. They exceed (in the NRS mode) T = lO"^ 
s in Cn and Z = 114 and T = lO" s in Z 120. 

In the neutron rich Hs isotopes the NRS barriers ex- 
tend up to relatively large values of the quadrupole mo- 
ment. Therefore their transition probabilities are smaller 
than in the RS mode. This fact together with the long 
a decay half-lifes implies that in this region one can find 
isotopes with very long half-lives. The RS fission is the 
dominant decay mode. However these very neutron rich 
Hs isotopes are extremely difficult to synthesize using 
contemporary experimental techniques. 



C. The region N > 180 

The saddle point of the first hump of the barrier takes 
its maximal energy value at A^ = 182. We can observe 
in these isotones another very high maximum of the fis- 
sion half-lives in all elements. In contrast, the A^ = 180 
isotones decay a few orders of magnitude faster than the 
neighboring nuclei with A^ = 178 and A^ = 182. The in- 
fiuence of the oblate minimum and the high first barrier 
is not strong enough to enlarge the fission half-lives here. 
In two isotones with A^ = 180, namely Cn and Z = 114, 
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the NRS fission is the dominant process with half-hves 
shorter than for the a decay. 

The energy difference between the oblate and the pro- 
late minima as well as between the absolute values of 
their quadrupole moments shrink in A'^ = 180 and N = 
182 continuously. Finally, one founds that all N = 184 
isotones do have a spherical ground state. This indicates 
a magic number at A'' = 184. Nevertheless, this feature 
does not have a big impact on the fission half-lives and 
the region of the most stable nuclei is slightly shifted 
towards the neutron deficient isotopes. 

At TV = 178 and iV = 182 the a decay half- lives also 
reach their maxima although they are less pronounced 
than for the fission half-lives. They correspond to the 
minima of Qa that can be observed for these neutron 
numbers, see Fig. 5. For the Z > 116 elements the a 
emission is the fastest decay process for isotopes lighter 
than N = 184. Most of them live long enough to allow 
the synthesis of these nuclei. The extremely long fission 
half-lives for the N = 178 and N = 182 isotopes arc 
larger than the a decay half-lives in Cn and Z = UA 
whereas for N = 180 NRS fission is the dominant mode. 

The isotopes with neutron number larger than the 
magic A'" = 184 are characterized by a rapid decrease 
of the fission half-lives with increasing neutron number. 
This is a consequence of the decrease of c^iwagy along the 
whole energy barrier. The second minimum goes down 
below the ground state and the height of the first barrier 
reduces substantially when heavier nuclei are considered. 
These strong trends can not be balanced by the few MeV 
decrease of the ground state energy due to the octupole 
deformation. Consequently, for the elements from Ds to 
Z = 120, fission becomes the dominant decay mode with 
half- lives below T = 10 /is. The a decay half-lives do 
not vary strongly along the isotopic chain in this region. 
The fluctuations are associated with the changes of the 
ground state deformations of the parent and the daugh- 
ter nuclei. In this region the decrease of the a decay 
half- lives with the proton number is also visible. Most 
isotopes in this region can not be synthesized due to the 
very short fission (in the proton deficient nuclei) or a 
emission (in the proton rich isotones) half-lives. The ex- 
perimental limit of T = 1 ms is exceeded only for a few 
nuclei. 

To conclude this section we would like to present in 
Fig. 15 the shortest half-lives of each isotope in the form 
of the chart of SHE's. From this figure it is easy to distin- 
guish the regions where each decay mode plays the most 
important role. The predominant decay mode, specially 
in the proton rich region, is the a emission. Roughly 
speaking, the spontaneous fission in the RS mode is dom- 
inant in the proton deficient nuclei with Z < lOA for 
N > 158, Z < 108 for N > 170, and Z < 116 for 
> 186. In the central part of the diagram the re- 
gion with the fastest decay in the NRS fission channel is 
deflned by 108 < Z < 114 and 170 < AT < 180. 

Two regions of long living nuclei can be found also in 
Fig. 15. The first one includes nuclei around ^^^Sgi62 



where half-lives reach T = 10** s. The other "island of 
stability" is associated with the anomalous long fission 
half-lives at N = 178 and A^ = 182. At these neutron 
numbers (and also at A^ = 184) long living isotopes of 
Cn and Z = 114 elements can be found. The half-lives 
of these isotopes are also longer than T = 10^' s. Very 
long half-lives characterize also the Hs isotopes decaying 
through the RS fission. 



VI. SUMMARY AND CONCLUSIONS 

In this work the Hatree-Fock-Bogoliubov theory with 
large basis size and the density dependent Gogny force 
as interaction have been used to study the most rele- 
vant properties of 160 heavy and super-heavy elements 
as well as their predominant decay modes. In order to 
keep the calculations as general as possible we allow for 
wave functions with different symmetries, namely axi- 
ally symmetric and reflexion symmetric, ax;ially symmet- 
ric and non reflexion symmetric and lastly triaxial wave 
functions. After a thorough analysis along different fis- 
sion paths and considering that our calciilations are ex- 
tremely CPU time demanding due to the large basis to 
study 160 nuclei we perform most calculations in the ax- 
ially symmetric approaches. 

The ground state deformations 132, Ps, (3i and p% as well 
as Qa factors, pairing properties and the two neutron and 
two proton separation energies are thoroughly discussed. 
The single particle energies are used as a guide for the 
interpretation of these properties. 

The fission paths for three representative nuclei, 
namely ^''^Hs, -^^^Ds and ^*^Cn, are analyzed in detail. 
Properties like mass parameters, pairing features and the 
variation of the action among others are calculated along 
the fission path with wave functions of the three types 
mentioned above. We find that though the shape of the 
fission barrier has a large impact on the fission half-lives, 
the mass parameter also plays an important role. In 
such a way that wave functions with larger barriers and 
smaller mass parameter may tunnel easier than alterna- 
tive ones with smaller barriers and larger mass. Since in 
general axially symmetric wave fimctions do have smaller 
masses than the triaxial ones, the restriction to axial sym- 
metry is a good option to perform a systematic study of 
the half-lives of SHE. In the second part of the paper a 
thorough study of the shapes of the barriers in the AS- 
RS and AS-NRS is performed. A reasonable explanation 
of heights and shapes as well as of octupole effects is 
obtained on the basis of the single particle energy lev- 
els. The two dimensional {Q2,Q?,) energy contour plots 
for the nuclei ^''**Hs and ^^^Cn allow to disentangle the 
different fission paths and its possible interconnections. 

In the third part of the paper the half-lives of all stud- 
ied nuclei are calculated for the different decay modes, 
namely a-decay and along the different fission paths. 
We find clear tendencies with the neutron number easily 
explainable on the basis of the behavior of Qa factors 
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and barrier shapes, respectively. In particular, we find 
that the a emission is the predominant decay mode, spe- 
cially in the proton rich region. Concerning the sponta- 
neous fission we obtain that the RS mode is dominant 
for the proton deficient nuclei in medium mass SHE's 
with Z < 104 for N > 158, Z < 108 for N > 170, and 
in the region of the heaviest SHE's with Z < 116 for 
A'' > 186. The fastest decay in the NRS fission channel 
takes place for 108 < Z < 114 and 170 < iV < 180. The 
long living nuclei can be found two regions. The first 
one is in the vicinity of ^^^Sgi62 where half-lives around 
T = 10^ s are obtained. The anomalous long fission half- 
lives for proton deficient isotones with N = 178, N = 182 
and N = 184 create the second region. At these neutron 
numbers several isotopes of the elements Cn and Z = lU 
are found with half-lives longer than T = 10^ s. Very 
long half-lives characterize also the Hs isotopes decaying 
through the RS fission. The nuclei beyond Z = 120 and 
A'' = 184 have half-lives too short to be detected within 
contemporary experimental limit of T = 10 /is. 

In conclusion we have presented a systematic study of 



SHE within the self-consistent HFB with very general 
wave functions and large configiiration space. Our calcu- 
lations provide an overall interpretation of the systcmat- 
ics and global properties of these elements and its decay 
modes. In general we find a reasonable good agreement 
with the experimental data. 
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